Electron-electron scattering and nonequilibrium noise in Sharvin contacts 
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We consider wide ballistic microcontacts with electron-electron scattering in the leads and calcu- 
late electric noise and nonlinear conductance in them. Due to a restricted geometry the collisions 
of electrons result in a shot noise even though they conserve the total momentum of electrons. We 
obtain the noise and the conductivity for arbitrary relations between voltage V and temperature T. 
The positive inelastic correction to the Sharvin conductance is proportional to T at low voltages 
eV <C T, and to \V\ at high voltages. At low voltages the noise is defined by the Nyquist relation 
and at high voltages the noise is related with the inelastic correction to the current by the Shottky 
formula Sin = 2e Un . 



Nonequilibrium electric noise is observed in most meso- 
scopic systems. It depends on the conduction mechanism 
and is more sensitive to the effects of electron-electron 
interactions than the average conductance In this 

article we are concerned with Sharvin-type ballistic con- 
tacts. In the absence of scattering near the contact, all 
the relaxation processes that lead to dissipation and a 
finite resistance of the contact take place deep in the 
leads, where the electron distribution is almost equilib- 
rium. As the motion of electrons in the nonequilibrium 
region near the contact is purely deterministic, the noise 
does not depend on the voltage and is specified by the 
Nyquist relation involving the equilibrium Sharvin con- 
ductance. If any impurities arc present in the contact, 
this results in a positive contribution to the resistance 
and in a shot noise, which is proportional to the cur- 
rent. Unlike the impurity scattering, electron-electron 
collisions do not contribute to the resistivity of a homoge- 
neous conductor with a parabolic spectrum because they 
conserve the total momentum of electrons. However very 
recently, it was shown both experimentally 0] and theo- 
retically Q that electron-electron scattering may result 
in a negative correction to the resistance of wide ballis- 
tic contacts. Therefore it is of interest to calculate the 
voltage-dependent electric noise in them and to find out 
whether collisions of electrons result in a shot noise like 
impurity scattering. 

Effects of electron-electron interaction on the shot 
noise have been extensively studied in the past for con- 
tacts with imperfect transmission. More than a decade 
ago, they were considered semiclassically for diffusive 
multichannel microbridges More recently, a num- 
ber of authors considered interaction effects in the shot 
noise of microstructures by modeling them as conduct- 
ing quantum dots that were either in the Kondo [5, 6] or 
Coulomb-blockade 0, S] regime. The interaction was as- 
sumed to take place between electrons in localized states 
on these dots. Naturally, this interaction strongly dif- 
fers from that in the bulk of the conductor. However 
our recent results show that even collisions of electrons 
far from the contact affect the average current and hence 
may cause its fluctuations. 

To calculate the noise, we use the semiclassical 



Boltzmann-Langevin method Previously, Kulik and 
Omclyanchuk used a similar approach to calculate elec- 
tric noise in Sharvin contacts caused by electron-phonon 
scattering in the zero-temperature limit [lo| . Here we 
extend this approach to arbitrary temperatures. 

We adopt a model of a ballistic contact similar to that 
of Kulik et al. [ll| for the case of electron-phonon scat- 



tering. Consider two 2D electron gases separated by a 
thin impenetrable barrier with a gap of width 2a. We 
assume that a is much larger than the Fermi wavelength 
and the screening radius but much smaller than both 
elastic and inelastic mean free path of electrons. The 
distribution functions of electrons on both sides of the 
insulator obey the Boltzmann equation 
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where E = — V(p is the electric field. The electron- 
electron collision integral in this equation is given by 



i [-/(P'k' ^ Pk) - J(pk ^ p'k') 



(2) 



p'kk' 



where 



J(pk ^ p'k') 



= iy(pk|p'k')/(p)/(k) 
x[l-/(p')] [l-/(k')] 



(3) 



and iy(pk|p'k') = Stt^ 



ek') '^(P + k — p' — k') is the probability of a transition 
from the state (p, k) to the state (p', k'). Here aee is a di- 
mensionless parameter of electron-electron scattering and 
1/ = m/iT is the two-dimensional density of states. Equa- 
tion dD) should be solved together with the Poisson equa- 
tion for the electric potential ip. It is possible to avoid 
solving the latter using the condition Ep 3> max(eV, T) 
[lot , which means that in the absence of collisions, elec- 
trons near the Fermi surface just move along straight 
lines. This condition allows us to set v = vpp/p and 
remove the term with electric field from Eq. ([T]). 

Now we should specify the boundary conditions to cal- 
culate the distribution functions using Eq. ([T]). We set 
/(p) = /o(ep) and ip ~ ±V/2 far from the gap in the left 
and right half-planes. 
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If collisions are neglected, /(p,r) depends solely on 
whether the electron trajectory originates from gap or 
not. It is convenient to use a notion of the angular do- 
main r2i„(r) that contains all the momenta of electrons 
that came to point r from the contact. In terms of this 
domain, the zero approximation distribution function is 

_ / /o(ep + e(^(r) T el//2), p i 17™(r) 

e(^(r)±ey/2), p G 17™(r) ^ ^ 

for the electrons in left (upper sign) and right (lower 
sign) half-spaces, respectively. 

The current through the contact is given by 



/ = e / dp 



(5) 



where v± is the component of v normal to the insula- 
tor and vector p ~ e||p labels points within the gap 
in the plane of insulator. Substituting expressions (|4]) 
into Eq. ([5]) results in the well known expression for the 
Sharvin conductance Gp = e^pFa/T^"^ ■ The first-order 
correction in scattering to Gq can be calculated by sub- 
stituting /(p, r) from Eq. ([T]) into Eq. ^ 



dp 



vj_ / dr/ee{/(°nP(r),r(T))}. (6) 



Here r is the time of travel to point p along the trajectory. 

The collision integral ([2]) involves four electron mo- 
menta p, k, p', and k'. If none of electrons with these 
momenta crosses the gap (i.e. falls within f2i„(r)), a sub- 
stitution of distribution functions ffl into ([51) results in 
lee = 0. As it was shown in Refs. [3| and |12| . the main 
contribution to the current ([6]) comes from collisions at 
points r located much farther from the gap than its size 
a. Hence ilm(r) may be considered as small and the 
fewer of the four momenta are in fli„ (r), the larger the 
contribution to the current. Also in Refs. 0] and [l2j 
it was shown that the largest contribution to the current 
comes from the collisions of electrons incident on the gap 
with electrons that are injected from the other half plane 
and have nearly opposite momentum. Therefore we can 
assume that only the electrons with momentum k are in- 
jected and lie in f2i„(r) while the electrons with the rest 
of momenta p, p', and k' are native to the considered 
half-plane. 

We sequentially integrate in Eq. © over the time, 
coordinate and momenta as it was done in Ref. 121 for 



the case of low voltage eV <C T. Here we consider the 
case of arbitrary voltages and obtain the correction to 
the current in a form of an integral over the energies 



ea'^aeem , Ic , , 



deic / dew / det.> 



2(2^)3 '"a 7 "^'^i 
X Fo(ep,ek,£p',£k') ^(£p + £k - £p' - £k') 

xe{D)/VD, (7) 



where Ic is a cutoff length much larger than 2a, which 
may be due to a very weak electron-impurity scattering 
or a finite size of the electrodes, 

Fo{ep, Ek, Ep', Ek') 

= [1 - hiSp)] [1 - /K(£k)] fUSp') /L(£k') 

- fUsp) fRie^) [1 - hisp')] [1 - /L(ek')], (8) 



and 



[(ep-ek)'-(£p'-£k')'] /4 



(9) 



is a value characterizing the deviation of the energies 
from the Fermi surface, which vanishes when all the en- 
ergies lie exactly at the Fermi surface. This expression 
allows us to analytically obtain the results in the limiting 
cases of high and low voltages and numerically calculate 
the correction for arbitrary relations eV/T. At high volt- 
ages eV/T 3> 1 the correction to the current has a form 



ea aeem , I 



(27r)3 



In - X (1 

a 



sign{V){eVf (10) 



and at low voltages eV/T <C 1 



/.n = ^^ln-x— (eT/)T (11) 

where the constant Gio = 3.72. 

Using the above semiclassical model, we can calculate 
the noise spectral density. It is expressed through the 
Fourier transform of the current correlation function as 
follows 



(12) 



We will calculate the spectral density at zero frequency 
1X1 = 0. Current fluctuation can be expressed in terms of 
fluctuation of the distribution function by the Eq. ([5]), 
so the current correlator has a form 



(Slit) 51(0)) = / dpi / dp 



(27r)2 J (27r)2 

— a —a 

X vij_V2i_ ((5/(pi,pi,t)(5/(p2,p2,0)). (13) 

The fluctuation df{p,r,t) obeys the Boltzmann-Lange- 
vin equation 



d d d \ df 

= (5/ee(p,r,t)+W^^*(p,r,i), 



(14) 



where i5J'^^*(p, r, <) is a Langevin source. The correlator 
of Langevin sources was calculated in |15| on the assump- 
tion that each collision is correlated only with itself and 
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equals 



(W^^*(ri, ti, pi) <5J^"*(r2, t2, P2)) - 2 "^(^i ~ ^2) 

X S{ti — (^pipa ^ (>/p'k'^pik + ^pik-»p'k') 
p'kk' 

+ ^ ^ ('/p'k'-i-pip2 + >^pip2-i-p'k') 
p'k' 

- 2^(Jpik^k'p2 + >/k'p2->pik) (15) 

kk' ^ 

We can neglect the field terms in for the same reason 
as it was done for the Boltzmann equation ([l}. As we 
consider finite temperatures, we have to take into account 
equilibrium fluctuations of / far from contact. To this 
end, we present 6f as a sum of fluctuation (5/o that has 
arrived from the depth of electrode by freely propagating 
without scattering [13[ and the integral of the right-hand 
side of Eq. p4|) over time t of travel to point r along the 
trajectory of a free electron 

00 

(^/IPiT,*) = (5/o(p,r,t) + J dr [(5/ee(p,r - vr, i - r) 



<5J"^*(p,r- vr,t-T)] . 



(16) 



Therefore the correlator of distribution functions in Eq. 
(fT3)) is given by 

(5/(pi,Pl,i)J/(p2,P2,0)) 

= (^/o(Pi,Pi,i)'^/o(P2,P2,0)) 

dr ((5/o(pi,Pl,i)(5/ee(p2,P2 - V2T, -r)) 
+ (<5/o(P2,P2,0)(5/ee(pi,pi - ViT, t - t)) 

00 00 

+ J '^'^l J dT2{5J^''\^l,pi--VxTl,t-Tl) 


X <5J'=-*(p2,P2-V2T2,-r2)). (17) 

Note that (S/o and 5 J*^^* are totally uncorrelated because 
of the causality principle. The first term in (17) is the 
two-time correlation function of fluctuations that origi- 
nate from the depth of electrodes and propagate to the 
point of observation without scattering. This correlation 
function is well known [l^ 

((5/o(pi, Pl, i)5/o(p2, P2, 0)) = (2^)2 5(pi - P2) 

x<5(pi - p2 - vit) /(pi)[l - /(pi)]. (18) 



Substituting this correlator to (|T3|) and (|T2|) results in the 
Nyquist equation Sa = 4rG'o, where Go is the Sharvin 
conductance. 



To the first order in the scattering, the spectral den- 
sity Sin ~ S — Sq is defined by the last three summands 
in Eq. (fT7|) . In the second and the third summands 
of this equation, 5Iee is obtained by variating the colli- 
sion integral lee @ with respect to 6f . To obtain the 
results to the first order in the interaction, we substi- 
tute the distribution functions ([4]) in 8Iee and evaluate 
the correlators {SfoSIee{Sf}) and {SIee{Sf} 5fo) using 
the zero-approximation correlators Eq. (fT8|) . The fourth 
summand in Eq. (fT7|) is obtained directly from Eq. (jlSp . 
Then we substitute the resulting expressions into Eqs. 
(HSl) and (HH). 

After some rearrangements the first-order spectral den- 
sity takes up a form 



00 a 



Sm = 2e Vyoi dt dpi / dp: 



d^pi f d'^p2 



(27r)2 J (27r)2 
X vi±V2± J dri J dT2 5{pi - P2 + V2T1 - vin - ^^2t) 



—(1 —a 
00 00 





6{t - n - T2) Ti + 5{t - n + T2) r2 



Pl-ViTi 



(19) 



where Fi and r2 are combinations of distribution func- 
tions and scattering fluxes 



1P2 



p'k' 

+ [1 - /(P2)] >/piP2^p'k' } + 2 X!{/(P2) •^Plk-.psk' 

kk' 

+ [1-/(P2)] Jp2k'^Pik}, (20) 

^2 = ^piP2[l - 2/(Pl)] ^ (-'^p'k'-^p2k - -'^pak^p'k') 
p'kk' 

+ [1 - /(Pl) - /(P2)] X!(Vk'^PlP2 - -^piP2^p'k') 

p'k' 

+ 2 [/(P2) - /(Pl)] ^(Jpik^p2k' - Jp2k'^pik). (21) 

kk' 

The two terms in Eq. (|19p have different physical mean- 
ing. The first of them corresponds to the case where the 
collision takes place during the time interval between the 
two crossings of the gap by the participating electrons 
(Fig. [T^). It originates from the collision integral in Eq. 
([T4l) and determines the Nyquist noise at ^ = 0. This 
term vanishes at T = for any V because it results from 
the equilibrium fluctuations in the depth of electrodes. 
The second term corresponds to the case where both 
electrons cross the gap after the collision (Fig. [TJd) and 
results from the corrections to the one-time correlation 
function of 5f. In equilibrium this function is a thermo- 
dynamic quantity and does not depend on the strength 
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FIG. 1: Fig. 1. Illustration of the two terms in Eq. ((T9j. An 
electron with momentum pi crosses the gap at time t > 0. 
Another electron with momentum p2 crosses the gap at time 
0. The collision takes place sd, t — ti < t. (a) The first term 
corresponds to the collision of electrons between the crossing 
of the gap; T2 is the time between the first crossing and the 
collision, (b) The second term corresponds to the collision 
before both crossings; T2 is the time between the collision and 
the first crossing. The inset shows the domains r2i„(r) and 
fiout(r). 



of scattering. Therefore the scattering corrections to it 
and the second term in Eq. p9|) vanish at y = 0. 

Consider now (|19p and isolate the dominant terms in 
it. To give a contribution to SFs, both the momenta pi 
and p2 must he either in the angular domain flin or in the 
centrally symmetric domain flout (see Fig. [TJ inset). As 
well as for the correction to the current, the contribution 
to the spectral density of noise is dominated by electron 
collisions far from the contact, so the angular domains are 
small and the contribution to the noise is maximum if a 
minimally possible number of electron momenta involved 
in scattering is restricted to or flout- On the other 
hand, the colliding electrons must have almost opposite 
momenta to ensure maximum phase space available for 
the scattering. 

In the first term with S{t ^ ti — T2), momenta pi and 
P2 have opposite directions (see Fig. [TJi) and lie in ilout 
and Qim respectively. To ensure maximum phase space 
for the scattering, pi and p2 must correspond either to 
the two initial or the two final states. Hence only the 
first sum should be retained in Eq. (pUl) . 

In the second term with (5(t — ri +T2), momenta pi and 
P2 are both in Qout (r) (see Fig. ^) and cannot have op- 
posite momenta. To make Eq. (PT|) nonzero, the electron 
with k must be in f2i„(r). Hence the dominant contribu- 
tion to Eq. (PT|) arises from the first term because Sp^p^ 
in it lifts one of the three restrictions on the momentum 
integration. 

We substitute the corresponding terms of Eqs. ([20)1 
and (pi]) into Eq. ([T5|. sequentially integrate over times, 
momenta and coordinates and obtain the spectral density 
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FIG. 2: Fig. 2. Dependencies of corrections to the current, 
measured in ea^aeem./(27r)^, and spectral density, measured 
in e^a^aeem/(27r)^ on eV/T. 



in a form of an integral over energies 



c^fl^Q; ui If f f f 
(2^)3 ^"f J J J dep, J de^, 

X 5{ep + Ek - Ep' - ek') Q{D) D-^'^ 
[i^i(ep,ek,£p',£k') + ^2(ep )], (22) 



where 



./fl(£k)[l - ./flXek)]{[l - fdep)] /L(ep')/L(ek') 
+ /L(ep)[l-/L(ep')][l-/L(ek')]} (23) 



and 



F2 = [l-2fL{ep)]Fo 



(24) 



with Fo defined by Eq. ©. 
At low voltages < T 



Si; 



(27r)3 



In^ X [2CioT2 + C2(eF)2] (25) 



where the constants Cio = 3.72 and C2 = 0.22 were 
calculated numerically. In view of Eq. ([TT]) , this is in full 
agreement with the Nyquist theorem. 

At high voltages eV ^ T the spectral density has a 
form 



(27r)3 



In-^ x21 (el-) 



(26) 



and is related with the inelastic contribution to the cur- 
rent (fTO]) by the Shottky formula Sin = '2elin- This is 
a consequence of the first approximation in the scatter- 
ing and the fact that the inelastic correction to the cur- 
rent is dominated by collisions far from the contact. The 
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weak scattering suggests that different electron collisions 
may be considered as independent random events whose 
contributions to the current simply sum up. As the colli- 
sions take place far from the contact, the angular domain 
riin is small and scattering of electrons within it may be 
disregarded. Therefore any scattering event changes the 
number of electrons crossing the contact by unity and 
this results in the classical shot noise of the inelastic cor- 
rection to the current. 

For arbitrary relations between voltage and tempera- 
ture the Sin{eV/T) and Iin{eV/T) dependencies can be 
obtained by numerically integrating Eqs. ([7]) and p2|. 
The resulting curves are shown in Fig. [2] 

To summarize, we have calculated the nonlinear correc- 
tion to the current and noise from electron-electron scat- 
tering for arbitrary relations between voltage V and tem- 
perature T. Both quantities are dominated by electron 
collisions at distances from the contact much larger than 
its size and are positive for all V. This is markedly dif- 
ferent from the case of impurity scattering, which results 
in a negative correction to the conductance and a cor- 
rection to the noise that is negative at low voltages and 
positive at high voltages. At low voltages, the correction 
to the noise is determined by thermal fluctuations that 
emerge from the depth of electrodes. At high voltages, 
it is determined by random collisions of nonequilibrium 
electrons and is related with the nonlinear correction to 
the current by the classical Shottky formula. 

An experimental test of the Shottky relation for wide 
ballistic contact in high-mobility samples could addition- 
ally verify that positive magnetoresistance and linearly 
increasing with temperature conductance observed in [2| 
are associated with electron-electron scattering at large 
distances from the contact. 
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